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(difference between two successive RR intervals) collected during atrial 
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estimated by the Kolmogorov-Smirnov test. If the lest density histogram 
is not significantly different from the standard density histogram, atrial 
fibrillation is delected. 
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1 

METHOD AND SYSTEM FOR DETECTION 
OF CARDIAC ARRHYTHMIA 

5 

BACKGROUND OF THE INVENTION 

1. Field of the invention: 

10 

The present invention relates to a method and a system for 
detecting cardiac arrhythmias from internally and/or externally detected 
activity of the heart . 

15 

2. Brief description of the prior art: 

Atrial fibrillation is a serious and common cardiac arrhythmia. Atrial 
fibrillation is associated with rapid, irregular atrial activation with life 

20 threatening sequelae such as stroke. The atrial activations are irregularly 
transmitted through the atrioventricular node leading to a correspondingly 
irregular sequence of ventricular activations as monitored by the 
ventricular interbeat (RR) intervals on the surface electrocardiogram 
(ECG). An RR interval is an interval between two successive heart beats. 

25 Clinically, in the surface ECG, atrial fibrillation is diagnosed by absence 
of P waves (normally associated with- the near synchronous activation of 
the atria) and a rapid irregular ventricular rate. P waves are difficult to 
determine automatically and irregular baseline activity of the ECG is 
common in atrial fibrillation. 

30 
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Although a number of different methods have been proposed to 
test for atrial fibrillation based on knowledge of the RR intervals and/or 
the surface ECG, the detection of atrial fibrillation based on this data ' 
nevertheless poses substantial problems (Murgatroyd, et al. "Identification 
5 of Atrial Fibrillation Episodes in Ambulatory Electrocardiographic 
Recordings: Validation of a Method for Obtaining Labeled R-R Interval 
Files," Pacing and Clinical Electrophysiology, (1995), pp.1 31 5-1 320). In 
the following description, the main strategies that have been proposed to 
assess atrial fibrillation based on knowledge of the RR intervals and/or 
10 ECG will be briefly reviewed. 

Since RR intervals during atrial fibrillation have a larger standard 
deviation and a more rapid decay of the autocorrelation function, there 
are proposals that the standard deviation and the autocorrelation function 

15 can be used to distinguish atrial fibrillation from sinus rhythm (Bootsma, 
et al. "Analysis of RR Intervals in Patients with Atrial Fibrillation at Rest 
and During Exercise," Circulation, (1970), pp.783-794). Since other 
abnormal rhythms also have a large standard deviation of RR intervals 
and a rapid decay of the autocorrelation function, these methods are 

20 difficult to apply in concrete situations. 

Moody and Mark (G. Moody, et al. "A New Method for Detecting 
Atrial Fibrillation Using R-R Intervals, 1 ' Computers in Cardiology, (1983), 
pp.227-230) classify RR intervals as short, long or regular. They then 

25 construct a Markov model for the probabilities for transitions between RR 
intervals in each of the three different length classes. Atrial fibrillation data 
has typical transition probabilities not shared by normal rhythms or other 
arrhythmia. Although the Markov model has high sensitivity for detecting 
atrial fibrillation, it tends to have a relatively low predictive value of a 

30 positive test. 
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Pinciroli and Castelli have investigated the morphology of 
histograms of RR intervals collected during atrial fibrillation and other 
arrhythmia (F. Pinciroli, et al. "Pre-clinical Experimentation of a 
5 Quantitative Synthesis of the Local Variability in the Original R-R Interval 
Sequence in the Presence of Arrhythmia," Automedica, (1986), vol.6, 
pp.295-317.Pinciroli and Castelli, 1986). They demonstrated that the 
histograms of the ratio between successive RR intervals show 
characteristic differences between normal rhythm and atria! fibrillation. 
10 The histogram of the ratio between successive RR intervals is 
symmetrical to the mean value. No quantitative methods were proposed 
to quantify the symmetry or to use it to develop a quantitative test. 

Since the baseline of the ECG is irregular during atrial fibrillation, 
15 Slocum (J. Slocum, et al. "Computer Detection of Atrial Fibrillation on the 
Surface Electrocardiogram," Computers in Cardiolody, (1987), pp. 253- 
254) has proposed that the regularity of the baseline, as determined by 
the power spectrum of the residual ECG after subtraction of the baseline 
of the QRS complexes can be used to detect atrial fibrillation. This 
20 method is necessarily sensitive to small amounts of noise that might 
corrupt the baseline of the ECG. 

Implantable ventricular and atrial defibrillators are devices that 
distinguish atrial and ventricular fibrillation from other rhythms. Typically, 

25 electrodes in these devices record intracardiac activity directly from the 
atria and ventricles. The methods that are used to detect atrial fibrillation 
in these devices cannot be easily applied to recordings that give 
information about the timing of the QRS complexes (U.S. Pat. No. 
6,144,878, issued to Schroeppel on November 7, 2000, U.S. Pat. No. 

30 6,035,233 issued to Schroeppel on March 7, 2000, U.S. Pat. No. 
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5,749,900 issued to Schroeppel on May 24, 1998, U.S. Pat. No. 
6,064,906 issued to Langberg et aL on May 16, 2000, U.S. Pat. No. 
5,772,604 issued to Langberg et al. on June 30, 1998, U.S. Pat.No. 
6,061,592 issued to Nigam on May 9, 2000, U.S. Pat. No. 5,951,592 
5 issued to Murphy on September 14, 1999, U.S. Pat. No. 5,941,831 issued 
to Turcot! on August 24, 1999, U.S. Pat. No. 5,591,215 issued to 
Greenhut et aL on January 7, 1997). 

Analysis of a histogram of the interbeat intervals can be used to 
10 discriminate between ventricular fibrillation and ventricular tachycardia. 
By counting the number of beats in predetermined interval classes, an 
algorithm identifies a given sequence as ventricular fibrillation or 
ventricular tachycardia (U.S. Pat. No. 5,330,508 issued to Gunderson on 
July 19, 1994). While this patent suggests that the invention is of value in 
15 detecting and treating atrial fibrillation (column 2, lines 29-31), it does not 
provide specific embodiment for detecting and treating atrial fibrillation. 

Based on the foregoing review of the prior art, it is apparent that 
there is a need to develop a method and a system for determining 

20 whether or not a given recording is atrial fibrillation based on the timing 
of the QRS complexes as measured from an internal and/or external 
monitor. Assessment of whether a patient is in atrial fibrillation based on 
the timing of the QRS complexes would be extremely useful, for example, 
for assessing the efficacy of specific drugs on a patient fitted with a 

25 monitoring device that measures the timing of the QRS complexes. 



30 
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5 

SUMMARY OF THE INVENTION 

In accordance with the present invention there is provided a 
5 method for detecting cardiac arrhythmia of a patient, comprising detecting 
RR intervals of the patient wherein each RR interval is an interval 
between two heart beats, constructing standard histograms of ARRs 
collected during cardiac arrhythmia of a plurality of subjects wherein each 
ARR is a difference between two successive RR intervals, constructing 
10 test histograms of ARRs of the patient from the detected RR intervals of 
. this patient, and comparing the standard and test histograms to detect 
whether the patient suffers from cardiac arrhythmia. 

In accordance with preferred embodiments: 

15 

- the standard and test histograms are probability density histograms; 

- a mean value of a given number of successive RR intervals of the 
patient is calculated, and a standard probability density histogram is 

20 chosen in relation to this mean value; 

-. the comparison of the standard and test histograms comprises 
adjusting a specificity-altering and sensitivity-altering parameter; 

25 - the comparison of the standard and test histograms comprises: 

- calculating a standard cumulative probability distribution 
from the standard ARR probability density histograms; 

- calculating a test cumulative probability distribution from 
the test ARR probability density histograms; 
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- computing a deviation between these standard and test 
distributions; and 

- detecting cardiac arrhythmia when the computed 
deviation is higher than the specificity-altering and 

5 sensitivity-altering parameter. 

The present invention also relates to a system for detecting 
cardiac arrhythmia of a patient, comprising: 

an RR interval detecting monitor detecting RR intervals of the 
10 patient, wherein each RR interval is an interval between two heart beat; 

a standard ARR histogram storage unit in which are stored 
standard histograms of ARRs collected during cardiac arrhythmia of a 
plurality of subjects, wherein each ARR is a difference between two 
successive RR intervals; 
15 a test ARR histogram calculator supplied with the detected RR 

intervals from the monitor and constructing test histograms of the ARRs 
of the patient; and 

a standard and test ARR histograms comparator supplied with 
the standard and test histograms, this comparator comprising a detector 
20 of cardiac arrhythmia of the patient responsive to the comparison of the 
standard and test histograms. 

It is within the scope of the present invention to apply the above 
concept to detection of not only atrial fibrillation but also to other cardiac 

25 arrhythmias including in particular but not exclusively atrial flutter, 
multifocal atrial tachycardia, ventricular tachycardia, premature ventricular 
contractions, etc., as well as to detection of other body phenomenon 
involving electrical activity. It is also within the scope of the present 
invention to use signals other than the RR intervals, histograms other 

30 than ARR probability density histograms, tests other than the KS test, and 
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7 

series of ARRs other than 100, and that other methods besides the 
Komogorov-Smirnov test can be used to compare test histograms with 
the standard histograms. 

5 The foregoing and other objects, advantages and features of the 

present invention will become more apparent upon reading of the 
following non restrictive description of a preferred embodiment thereof, 
given by way of example only with reference to the accompanying 
drawings. 

10 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 In the appended drawings: 

Figure 1 are time series showing the RR intervals from subject 202 . 
from the MIT-BIH arrhythmia database. The solid line directly under the 
time series of RR intervals shows the assessment of atrial fibrillation 
20 (indicated by AF) or non-atrial fibrillation (indicated by N) as reported in 
the database. The solid line at the bottom of Figure 1 indicates the 
assessment of atrial fibrillation, indicated by 1, and non-atrial fibrillation, 
indicated by 0, based on an algorithm presented herein. 

25 Figure 2 is a flow chart illustrating a preferred embodiment of the 

method according to the present invention, for detecting atrial fibrillation 
based on RR intervals. 

Figure 3a-3p are ARR standard probability density histograms 
30 during atrial fibrillation. Mean RR intervals are a) 350-399 ms, b) 400-449 
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ms, c) 450-499 ms, d) 500-549 ms, e) 550-599 ms, f) 600-649 ms, g) 
650-699 ms, h) 700-749 ms, i) 750-799 ms, j) 800-849 ms, k) 850-899 
ms, I) 900-949 ms, m) 950-999 ms, n) 1000-1049 ms, o) 1050-1099 ms, 
and p) 1100-1049 ms. 

5 

Figure 4 is the standard deviation of ARR which consists of the 
standard ARR probability density histogram as a function of mean RR 
interval. 

10 Figure 5 illustrates the Kolmogorov-Smirnov (KS) test. A 

cumulative probability distribution based on patient test data is compared 
with a standard cumulative probability distribution. D is the greatest 
distance between two cumulative probability distributions. 

15 Figures 6a and 6b show an example of the standard deviation 

(Figure 6a) and the skewness (Figure 6b) of a test ARR probability 
density histogram. The line represents the standard deviation (Figure 6a) 
and the skewness (Figure 6b) of the standard ARR probability density 
histogram. 

20 

Figure 7 shows the receiver operating characteristic curve (ROC) 
when this method is tested on the MIT-BIH atrial fibrillation/flutter 
database. The specificity increases with increase in P ct while the 
sensitivity decreases with an increase in P c . 

25 . 

Figure 8 is a block diagram of a preferred embodiment of the 
system according to the present invention for implementing the method 
of Figure 2, for detecting atrial fibrillation based on RR intervals. 
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Figure 9 is a block diagram of a preferred embodiment of a test 
and standard ARR histogram comparator forming part of the system of 
Figure 8. 

5 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Although the preferred embodiment of the present invention will be 
10 described in relation to atrial fibrillation, the same concept can be applied 
to detection of other cardiac arrhythmias including in particular but not 
exclusively atrial flutter, multifocal atrial tachycardia, ventricular 
tachycardia, premature ventricular contractions, etc. This concept can 
also be applied to detection of other body phenomenon involving 
15 electrical activity. 

Data was obtained from the MIT-BIH atrial fibrillation/flutter 
database. The data contains 300 atrial fibrillation episodes, sampled at 
250 Hz for TO hours from Holter tapes of 25 subjects. Arrhythmia 
20 detection was carried out by trained observers and was confirmed by an 
independent evaluation. 

Figure 1 is a typical time series of RR intervals from a patient with 
atrial fibrillation. Immediately under the recording is a solid marker line 

25 101. When atrial fibrillation occurs this marker line 101 is set to AF; 
otherwise it is set to N, which indicates a rhythm that is not atrial 
fibrillation. The graph of Figure 1 also shows a lower solid line 102 
indicating the assessment of atrial fibrilation, indicated by 1, and non- 
atrial fibrilation, indicated by 0, based on an algorithm according to the 

30 present invention. At the onset of atrial fibrillation, the rhythm dramatically 
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changes to irregular with large fluctuation. In paroxysmal atrial fibrillation 
there is sudden starting and stopping of atrial fibrillation. 

Figure 2 shows a flow chart of a preferred embodiment of the 
5 method according to the invention for detecting atrial fibrillation. Figure 8 
is a block diagram of a preferred embodiment of the system according to 
the invention for implementing this method. 

The standard ARR probability density histograms are prepared as 
10 described hereinafter before the detection of atrial fibrillation, and then 
stored in an adequate storage unit 804 (Figure 8). 

RR intervals of the patient are first detected (201 of Figure 2) 
through an internal and/or external RR interval monitor 801 (Figure 8) 
15 detecting electrical activity of the heart beat of the patient. 

ARR is defined as the difference between two successive RR 
intervals. In the preferred embodiment, blocks of 100 successive RR 
intervals are processed during atrial fibrillation. For that purpose, the 
20 detected RR intervals from the monitor 801 are counted (202 of Figure 2) 
by a RR interval counter 802 (Figure 8) until the number of detected RR 
intervals reaches 100 intervals (203 of Figure 2). 

The mean value of each block of 100 RR intervals is computed 
25 (204 of Figure 2) by means of a calculator 803 from the RR intervals from 
the monitor 801. Of course, the calculator 803 is supplied with the count 
from the counter 802. This mean value identifies the block of 100 RR 
intervals as falling into one of sixteen (16) different classes, respectively 
corresponding to mean values of RR between 350-399 ms, 400-449 ms, 
30 450-499 ms, 500-549 ms, 550-599 ms, 600-649 ms, 650-699 ms, 700- 
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749 ms, 750-799 ms, 800-849 ms, 850-899 ms, 900-949 ms, 950-999 
ms, 1000-1049 ms, 1050-1099 ms, and 1100-1049 ms. For each of the 
sixteen (16) classes, a standard ARR probability density histogram is 
compiled by lumping data together from all the subjects, for example the 
5 subjects of the above mentioned MIT-BIH atrial fibrilation/flutter database. 
The resulting histograms (see for example in Figures 3a-3p) are taken to 
be the standard ARR probability density histograms for atrial fibrillation, 
sorted by the mean RR interval (see for example in Figures 3a-3p) and 
stored in storage unit 804. In other words, a standard ARR histogram 
10 selector 805 chooses the standard ARR probability density histogram 
(Figures 3a-3p) corresponding to the class in which the computed mean 
value of RR intervals (from 204 in Figure 2) of the block of 100 RR 
intervals under consideration falls (205 of Figure 2). 

15 Obviously, it is within the scope of the present invention to 

construct the standard ARR probability density histograms using a 
different number of consecutive RR intervals, for example 25, 50 or any 
other number of consecutive RR intervals. It is also within the scope of 
the present invention to construct the standard ARR probability density 

20 histograms using mean RR intervals that lie in other ranges, for example 
300-399 ms, 400-499 ms, 500-599 ms, etc. 

Figure 4 shows the standard deviation (SD) of the standard 
probability density histograms of ARR. 

25 

Test ARR probability density histograms are constructed (206 of 
Figure 2) by a calculator 806 from the data obtained from the patient (test 
record) through the monitor 801. As indicated in the foregoing description, 
the blocks of 100 successive RR intervals are determined by the counter 
30 802. In order to test for atrial fibrillation in a test record, the test ARR 
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probability density histograms based on the blocks of 1 00 successive RR 
intervals from the test record, are compared (207 and 208) through a 
comparator 807 to the chosen standard ARR probability density 
histograms from selector 805. 

5 

In the test ARR histogram calculator 806 a sequence of 100 RR 
intervals is centered on each beat in turn, and the relevant test ARR 
probability density histograms are calculated. A/so, a standard 
cumulative probability distribution is calculated by integrating the area 
10 under the curves of the standard ARR probability density histograms, and 
a test cumulative probability distribution is computed by integrating the 
area under the curves of the test ARR probability density histograms 
(Figure 5). 

15 The similarities between the test histograms for a given patient 

and the standard histograms are evaluated in the test and standard ARR 
histogram comparator 807 using the above mentioned Kolmogorov- 
Smirnov (KS) test (207 and 208 of Figure 2). As indicated, Figure 5 
shows an example of cumulative probability distributions of standard 

20 histograms (standard curve) and test histograms (test curve). 

Referring to Figure 9, a calculator 901 (Figure 9) computes the 
cumulative probability distribution of the standard probability density ARR 
histograms. A calculator 902 computes the cumulative probability 
25 distribution of test probability density ARR histograms. According to the 
KS test, one assesses if two given distributions are different from each 
other. In other words, the greatest vertical distance D between the two 
cumulative probability distributions is measured by a calculator 903 which 
returns a p value in the following manner: 
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where X = +0A2 + 0A1/ ffi)* D. n ^ ^~J^ • ^ is the 

number of data points on the standard cumulative probability distribution. 
5 A/ 2 is the number of data points in the test cumulative probability 
distribution. A detector 904 determines whether the p value is greater 
than a certain, appropriately selected threshold P c , and detection of p > 
P c indicates that the cumulative probability distributions are not 
significantly different from one another. Since the standard ARR 
10 probability density histograms is representive of atrial fibrillation, a value 
of p > P c constitutes a positive identification of atrial fibrillation (or more 
accurately failure to reject the hypothesis that the test cumulative 
probability distribution is not atrial fibrillation) (208 in Figure 2). 

15 Figures 6a and 6b show a comparison of the ARR probability 

density histograms in terms of standard deviation and skewness. A small 
D defined above indicates that the standard deviation and the skewness 
of a test ARR probability density histogram are clustered around those of 
the standard ARR probability density histograms. 

20 

The results were assessed by four categories as followed: true 
positive (TP) — atrial fibrillation is classified as atrial fibrillation; true 
negative (TN) — non-atrial fibrillation is classified as non-atrial fibrillation; 
false negative (FN) - atrial fibrillation is classified as non-atrial fibrillation; 
25 false positive (FP) - non-atrial fibrillation is classified as atrial fibrillation. 
Sensitivity and specificity are defined by TP/(TP+FN) and TN/(TN+FP), 
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respectively. The predictive value of a positive test (PV+) and the 
predictive value of a negative test (PV-) are defined by TP/(TP+FP) and 
TN/(TN+FN), respectively. 

5 The receiver operating characteristic curve (ROC) gives the 

sensitivity and the specificity* in the artrial fibrillation detection algorithm. 
Variation of the value of P c determines the ROC. Figure 7 shows the 
ROC of the assessment of the KS test for the MIT-BIH atrial 
fibrillation/flutter database. Reducing P C1 the sensitivity increases and the 
10 specificity decreases. Assuming P c = 0.003944, the sensitivity is 96.5%, 
the specificity is 96.5%, the PV+ is 95.2% and PV- is 97.5%. P c is 
therefore a sensitivity-altering and specificity-altering parameter. 

It will appear to those of ordinary skill in the art that the method of 
15 Figure 2 and the system of Figure 8 can be implemented through a 
properly programmed computer. 

Although the present invention has been described hereinabove 
by way of a preferred embodiment thereof, this embodiment can be 
20 modified at will, within the scope of the appended claims, without 
departing from the spirit and nature of the subject invention. 
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WHAT IS CLAIMED IS: 



1. A method for detecting cardiac arrhythmia of a patient, 
5 comprising: 

detecting RR intervals of the patient, wherein each RR interval is 
an interval between two heart beats; 

constructing standard histograms of ARRs collected during cardiac 
arrhythmia of a plurality of subjects, wherein each ARR is a difference 
10 between two successive RR intervals; 

constructing test histograms of ARRs of the patient from the 
detected RR intervals of said patient; and 

comparing said standard and test histograms to detect whether 
said patient suffers from cardiac arrhythmia. 

15 

2. A method for detecting cardiac arrhythmia as defined in claim 
1, wherein the standard and test histograms are probability density 
histograms. 

20 3. A method for detecting cardiac arrhythmia as defined in claim 

1, further comprising calculating a mean value of a given number of 
successive RR intervals of the patient, and choosing one of the standard 
probability density histograms in relation to said mean value. 

25 4. A method for detecting cardiac arrhythmia as defined in claim 

3, wherein choosing one of the standard probability density histograms 
comprises selecting the standard histogram corresponding to the range 
of RR intervals in which said mean value is located. 
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5. A method for detecting cardiac arrhythmia as defined in claim 
1 , wherein said comparing of the standard and test histograms comprises 
adjusting a specificity-altering and sensitivity-altering parameter. 

5 6. A method for detecting cardiac arrhythmia as defined in claim 

1, wherein said comparing of the standard and test histograms comprises 
calculating a standard cumulative probability distribution from said 
standard ARR probability density histograms, calculating a test 
cumulative probability distribution from said test ARR probability density 
10 histograms, and computing a deviation between said standard and test 
distributions. 

7. A method for detecting cardiac arrhythmia as defined in claim 
5, wherein said comparing of the standard and test histograms 
15 comprises: 

calculating a standard cumulative probability distribution from said 
standard ARR probability density histograms; 

calculating a test cumulative probability distribution from said test 
ARR probability density histograms; 
20 computing a deviation between said standard and test 

distributions; and 

detecting cardiac arrhythmia when the computed deviation is 
higher than said specificity-altering and sensitivity-altering parameter. 

25 8. A system for detecting cardiac arrhythmia of a patient, 

comprising: 

a RR interval detecting monitor detecting RR intervals of the 
patient, wherein each RR interval is an interval between two heart beats; 
a standard ARR histogram storage unit in which are stored 
30 standard histograms of ARRs collected during cardiac arrhythmia of a 
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#■ 

plurality of subjects, wherein each ARR is a difference between two 
successive RR intervals; 

a test ARR histogram calculator supplied with the detected RR 
intervals from the monitor and constructing test histograms of said ARRs 
5 of said patient; and 

a standard and test ARR histograms comparator supplied with said 
standard and test histograms, said comparator comprising a detector of 
cardiac arrhythmia of the patient responsive to the comparison of said 
standard and test histograms. 

10 

9. A system for detecting cardiac arrhythmia as defined in claim 8, 
wherein the standard and test histograms are probability density 
histograms. 

15 10. A system for detecting cardiac arrhythmia as defined in claim 

8, further comprising a calculator of a mean value of a given number of 
successive RR intervals of the patient, and a selector of one of the 
standard probability density histograms in relation to said mean value. 

20 1 1 . A system for detecting cardiac arrhythmia as defined in claim 

10, wherein the selector of one of the standard probability density 
histograms comprises means for selecting the standard histogram 
corresponding to the range of RR intervals in which said mean value is 
located. 

25 

12. A system for detecting cardiac arrhythmia as defined in claim 
8, wherein said standard and test ARR histograms comparator comprises 
an adjustable, specificity-altering and sensitivity-altering parameter. 
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18 

1 3. A system for detecting cardiac arrhythmia as defined in claim 
8, wherein said standard and test ARR histograms comparator comprises 
a calculator of a standard cumulative probability distribution from said 
standard ARR probability density histograms, a calculator of a test 

5 cumulative probability distribution from said test ARR probability density 
histograms, and a calculator of a deviation between said standard and 
test distributions. 

14. A method for detecting cardiac arrhythmia as defined in claim 
10 12, wherein the standard and test ARR histograms comparator 

comprises: 

a calculator of a standard cumulative probability distribution from 
said standard ARR probability density histograms; 

a calculator of a test cumulative probability distribution from said 
15 test ARR probability density histograms; 

a calculator of a deviation between said standard and test 
distributions; and 

a detector of a cardiac arrhythmia when the computed deviation is 
higher than said specificity-altering and sensitivity-altering parameter. 

20 
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